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Motivated by experimental limitations commonly met in 
the design of solid state quantum computers, we study the 
problems of non-local Hamiltonian simulation and non-local 
gate synthesis when only homogeneous local unitaries are per- 
formed in order to tailor the available interaction. Homoge- 
neous (i.e. identical for all subsystems) local manipulation 
implies a more refined classification of interaction Hamiltoni- 
ans than the inhomogeneous case, as well as the loss of univer- 
sality in Hamiltonian simulation. For the case of symmetric 
two-qubit interactions, we provide time-optimal protocols for 
both Hamiltonian simulation and gate synthesis. 

03.67.-a, 03.67.Lx 



I. INTRODUCTION 

A. Simulation of non-local Hamiltonians and 
synthesis of non-local gates 

The simulation of Hamiltonians has attracted consi- 
derable attention in several areas of quantum physics, 
including those of nuclear magnetic resonance JL | (NMR) , 
quantum control [Q and quantum information [ 3|-0 . Let 
us consider a quantum system that naturally evolves ac- 
cording to some Hamiltonian H for a time t. By inter- 
spersing its evolution exp(-iHt) with a sequence of fast 
unitary operations Wk B, another evolution 
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is obtained which, when \Htk\ << 1, approximately cor- 
responds to having the Hamiltonian 
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Vfc = Wk ■ ■ ■ Wi [we assume V n — I], acting for a time 
t' =t/c, where c > is the time overhead of the simula- 
tion. This means that we can effectively generate, given 
an available Hamiltonian H and some control operations 
Wk, a new Hamiltonian H' , which may otherwise not 
be accessible experimentally. Hamiltonian simulation is 
therefore a powerful technique to control the evolution 
of a quantum system. As such, it has a broad range of 



potential applications, and may eventually become com- 
monplace in experimental implementations of quantum 
information and quantum computation. 

The synthesis of a gate U = exp(-iH'T) by using H 
and control operations Wk is a closely related issue, for 
U can be achieved by simulating H 1 for a time T. We 
notice, however, that inequivalent sequences of control 
operations Wk may also lead to the same gate U — the 
only requirement is that at the end of the protocol the 
gate has been performed — and thus simulating H' may 
not be the most convenient way to proceed. 

A series of recent contributions l^-Q has recently ad- 
dressed, in the context of multipartite quantum systems, 
the problem of simulating one non-local Hamiltonian (or 
interaction) by using another non-local Hamiltonian and 
local unitary transformations (LU) as control operations. 
In particular, special emphasis has been placed on assess- 
ing the universality of this model , as well as on de- 
termining time-optimal simulation protocols and devel- 
oping a quantitative theory of interactions based on their 
simulation capabilities |^,|6| . Similarly, time-optimal pro- 
tocols for non-local gate synthesis have been thoroughly 
analyzed in the context of quantum control ^] , where the 
problem has been reduced, in the case of a two-qubit sys- 
tem, to a problem of non-local Hamiltonian simulation. 
In jjlj the corresponding optimal non-local Hamiltonian, 
and thereby the minimal time for two-qubit gate synthe- 
sis, have been presented. 



B. Hamiltonian simulation in solid state physics: 
homogeneous local manipulation 

In most of previous work on non-local Hamiltonian 
simulation, arbitrary LU are assumed to be available 
as control operations. This model is motivated by the 
fact that in a number of experimental schemes for quan- 
tum information processing (using e.g. ion traps, neutral 
atoms or photons |Tl|] ) independent LU operations can be 
performed on each of the interacting systems supporting 
the qubits. In solid state implementations, however, the 
relevant systems are frequently too close to each other to 
be addressed individually. A good example of this is a 
set of spins in a NMR molecule or in a lattice of quantum 
dots, where an extraordinarily focused (and thus unfea- 
sible) magnetic field would be required to independently 
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address each spin. Therefore, and due to experimental 
limitations, the current model for interaction simulation 
does not apply to a number of solid state schemes. 

Most relatedly, an alternative, solid-state oriented 
model for quantum computation has recently been put 
forward @. It overcomes the above lack of controllabil- 
ity by encoding logical qubits in convenient subspaces of 
few physical qubits, in a way that universal computation 
can be accomplished by just tunning the intensity of a 
specific non-local Hamiltonian such as the isotropic ex- 
change interaction, naturally available in many systems. 
Notice, however, that in solid state not all forms of local 
manipulation are equally inaccessible. In particular, a 
uniform magnetic field can be easily applied on all spins, 
inducing the same local evolution on each of them and 
thereby producing what we shall call an homogeneous 
local unitary transformation (HLU). In this sense, HLU 
appear as a rather restricted, but definitely feasible, set 
of control operations in solid state physics, that may be 
worth considering. 

C. Results 



A remarkable feature of the HLU model of simulation 
will be that it is ruled by the usual majorization rela- 
tion [Q, instead of the special majorization relation that 
rules the non-homogeneous LU model @. As we shall ex- 
plain, this gives a neat picture of which two-qubit Hamil- 
tonian simulations are possible in the symmetric case: 
H can simulate H' if and only if this implies a gain of 
isotropy in the interaction. 

The organization of this paper goes as follows. In sect. 
II we address the problem of two-qubit Hamiltonian sim- 
ulation and present the optimal solution in the case of 
symmetric and antisymmetric interactions. Using these 
results, the time-optimal synthesis of symmetric gates 
using a symmetric interaction is analized in sect. III. In 
sect. IV we discuss the synthesis of symmetric gates us- 
ing an arbitrary interaction. In sect. V we also show 
how to achieve any two-qubit gate by allowing for two 
non-homogeneous interventions. The paper finishes with 
some conclusions. 

II. OPTIMAL HAMILTONIAN SIMULATION 
USING HLU 



In this paper we address the problems of non-local 
Hamiltonian simulation and gate synthesis using a given 
non-local Hamiltonian and fast HLU. Our results, mostly 
concerned with two-body interactions in multi-qubit sys- 
tems, imply that an interaction, when only supplemented 
with HLU, does not constitute a universal resource for 
simulating other interactions or synthesizing arbitrary 
gates, in sharp contrast with the inhomogeneous LU 
model jl|,|]|7]]- Nevertheless, the HLU model still allows 
for a rich variety of Hamiltonian simulations and gate 
synthesis. In particular, HLU can be used to eliminate 
unwished anisotropic terms (as coming e.g. from a spin- 
orbit interaction) in a system of spins with imperfect 
isotropic exchange interaction. This task is of practical 
relevance, for instance, in the context of encoded uni- 
versality, and has been considered in jl3| . Therefore the 
HLU model for non-local Hamiltonian simulation pro- 
vides us with a complementary or alternative technique 
to that discussed in Jl3j . 

We shall present the following results. 

(i) Necessary and sufficient conditions for two-qubit 
Hamiltonian simulation under HLU, both for symmetric 
and antisymmetric interactions. These results also apply 
to some multiple-qubit scenarios with two-body interac- 
tions. 

(ii) Time-optimal synthesis of two-qubit gates using a 
symmetric interaction. 

(iii) Recovery of universality in gate synthesis using HLU 
and just one initial and one final non-homogeneous LU 
per two-qubit gate. 



Consider a collection of N interacting systems, say N 
qubits for concreteness, that evolve according to a Hamil- 
tonian H\2...n- If only fast HLU can be performed on the 
qubits in order to control their evolution, i.e. Wk = wf N 
in Eq. (|l|), then, according to Eq. (0), any simulated 
Hamiltonian H' 12 ...n must be of the form 

H' 12 ... N = cJ2pk(vf N )H 12 ... N (vf N )\ (3) 

k 

where pk = tk/t sum up to one. It already follows from 
this expression that, in general, it is not possible to sim- 
ulate an arbitrary Hamiltonian H' 12 ,.. N by H12...N and 
HLU alone. As an example, suppose that Hn- N is sym- 
metric under exchange of qubits. Then H' 12 ^ N neces- 
sarily has the same symmetry JT^|, so that we cannot 
simulate any inhomogeneous evolution. Thus, the HLU 
model for simulation is rather restricted. Nevertheless, 
it is still a useful and — most important — available tool 
to control the evolution of the qubits. 

We start our analysis by considering N — 2 qubits, for 
which we present most of our results. Later in this section 
we shall present some generalizations to the general N 
case. 

A. Two-qubit Hamiltonians 

Consider two qubits, A and B, that naturally evolve 
according to some Hamiltonian Hab ■ The Hamiltonians 
achievable by interspersing the evolution of Hab with 
control operations Wk <£> Wfe, Wk G SU(2), are now 
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h 'ab = c ^jPk v k ® v k H AB v\ ®v\, c> 0, (4) 

k 

where we have followed the same steps than in ||. Be- 
cause the set {/, <Ji, a 2, ct 3 } forms a basis for the self- 
adjoint operators acting on C 2 , we can always expand 
Hab as 

3 

H AB = al ® I + h A ® I + I ® h B + ( M ) 

ij <J -i ® (Tj , 

(5) 

where M is a real 3x3 matrix called the "Pauli repre- 
sentation" of Hab, and h A and h B are some local terms, 
that we set to h A — h B = h for simplicity []l6| . The term 
al ® / is an irrelevant phase and we ignore it. We can 
always add or eliminate equal local terms by applying an 
appropriate HLU, e~ lht ® e~ lht , because 

(e- lht ® e~ iht ) e- lHt = e -HMi+iah+B)t + Q ^ (6) 

Therefore, we can cancel the local terms of H AB , simulate 
the interaction part of H' AB using the interaction part of 
Hab and, then, add the local terms of H' AB with the cor- 
responding HLU. Hence, we need consider Hamiltonians 
with only interaction terms. For any unitary u € SU(2) 
there always exists a rotation R £ 50(3) such that 

uoiU* = RfjCj- (7) 
j 

It is worth mentioning that this argument only works for 
qubits |^7|| . Using this, we write (Q) as 

M' = cY d PkRkMRl . (8) 

k 

The resulting M' is thus proportional to an orthogonal 
(real unitary |L8|]) mixing of cM. Every term R^MR^ 
in Eq. (||) conserves the trace and the symmetry of M. 
This implies that 

tr(M') = ctr(M). (9) 

Thus, we have no freedom in choosing the time overhead c 
(inverse of the efficiency s of ||) unless tr(M)—tr(M') = 
0. Another consequence of Eq. (||), is that with a sym- 
metric (antisymmetric) M, M = M T (M = -M T ), we 
can only simulate a symmetric (antisymmetric) M' . Let 
us study separately these two cases. 

1. Symmetric M 

Let us start by considering a Hamiltonian symmetric 
under the exchange of the subsystems, SHabS^ = Hab, 
where the swap operator S is defined as to exchange the 



states of qubit A and B, S(\i(j)a <8 \4>)b) = \4>)a ® W)b- 
We note that this is the case of many relevant interac- 
tions, such as the Ising interaction a x ®a x , the XY model 
or anisotropic exchange interaction a x ® a x + a y ® a yi 
and the Heisenberg or isotropic exchange interaction 

<J X ® CF X + CTy ® CFy + (J z ® a z . 

Symmetric real matrices are the only ones that can be 
diagonalized with SO(3) transformations. Calling A (A') 
the vector of eigenvalues of M (M 1 ) we can rewrite (0) 
as 

k 3 

where we have absorbed the orthogonal matrices that 
diagonalize M and M' in redefinitions of the Rk matrices. 
Doubly stochastic matrices form a convex set with the 
permutation matrices as the generators (see Ref. |fL8|). 
The squared elements of an orthogonal matrix, (R)fj, 
always form a doubly stochastic matrix. Thus, when Eq. 
( |Io| ) is satisfied there exists a doubly stochastic matrix, 
D, such that 

A' = c DX. (11) 

On the other hand, for every permutation matrix P 
there exist orthogonal matrices R fulfilling (R)fj = (P)ij- 
Therefore, given a doubly stochastic matrix D we al- 
ways can find an orthogonal mixing such that (D)ij — 
J2k P k (Rk)ij- Concluding, M simulates M' iff there ex- 
ists one doubly stochastic matrix D fulfilling Eq. dO). 
The last is equivalent to say that cA majorizes A' 
(the proof is in jl8| ) . 

Result 1: Let Hab and H' AB be two symmetric in- 
teractions. Hab simulates H' AB by using fast HLU if, 
and only if, the eigenvalues of M' are majorized by the 
eigenvalues of cM , with the appropriate time overhead c 
given by Eq. ^). 

Using the notation of Q , we write Result 1 as 

H' AB <hlu cHab <=> A' -< cA. (12) 

Notice that the above simulation strategy is automati- 
cally optimal because the time overhead c is fixed. 

The above result can be understood in physical terms 
as follows. The more isotropic a Hamiltonian is, the less 
simulating power it carries. External homogeneous ma- 
nipulations can only produce mixing, thus isotropy. More 
precisely, the isotropy of Hab — as measured by the de- 
gree of mixing of A — sets the arrow of allowed simu- 
lations by HLU. Accordingly, highly anisotropic Hamil- 
tonians, e.g. the Ising model type Hab — &z ® o~ z , 
are excellent resources to simulate other Hamiltonians 
and, later on, to synthesize efficiently all sort of non- 
local gates. Instead, the isotropic exchange interaction 
Hab — o~ x ® <j x + o-y ® a y + a z ® a z is a fixed point under 
HLU transformations, and can not simulate any other 
interaction. 
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2. Time-optimal simulation with antisymmetric M 

We consider now the case of an interaction fulfilling 
SHabS^ = —Hab- There is an isomorphism between 
3x3 antisymmetric matrices and 3-dimensional vectors 
via v{M)i = eijkMjk- Due to the invariance of the Levi- 
Civita tensor e^fc under rotations, we obtain v(RMR T ) = 
Rv(M). Then, we can write Eq. (JsJ) in vectorial form: 

v(M') =cJ2Pk Rk v(M). (13) 

k 

By the triangular inequality we have the following rela- 
tion among the two moduli 

v(M') < c v(M). (14) 

Result 2: Let Hab and H' AB be two antisymmetric 
interactions. Hab can always simulate H\ B using fast 
HLU for any time overhead c fulfilling \14), i.e., 

H'ab <hlu cH A b v' <cv. (15) 



convenient HLU, we can make the system evolve accord- 
ing to first-neighbour anisotropic (or isotropic) exchange 
interactions. More generally, we can achieve any first- 
neighbour interaction compatible with Eq. (|T^). 

III. TIME-OPTIMAL SYNTHESIS OF 
NON-LOCAL GATES USING HLU 

The problem of producing a two-qubit gate U in opti- 
mal time using a given interaction Hab and fast LU was 
throughoutly explored in |J and finally solved in Jl(| . In 
this section we study this problem in the case where the 
given interaction Hab is symmetric (SHabS^ = Hab) 
and only HLU are allowed. Notice that symmetric inter- 
actions are the most common ones in nature. Because 
HLU are also symmetric, the net effect of interspersing 
the action of Hab with HLU must give a symmetric U. 
Thus, we only consider unitary gates U fulfilling 

SUS^ = U. (17) 



The minimal time overhead or optimal simulation is 
achieved when the equality ( |l4| ) holds, c = v' /v. An in- 
teresting consequence of Result 2 is that any Hamiltonian 
with antisymmetric M can simulate itself backwards in 
time with overhead 1, also that it can simulate the zero 
Hamiltonian. In other words, we can always stop its ac- 
tion. We shall use this property later on. 

We finish our analysis for two qubits by noticing that 
the above cases correspond to the general decomposition 
of a rank two tensor in irreducible representations of the 
rotation group, J = 0, J = 1 and J = 2. The scalar 
pieces produces the Eq. (^) restriction, the antisymmet- 
ric case comes from J — 1 and gives Eq. ([l3| ) and the 
standard majorization is related to the J = 2 irreducible 
representation. 

B. Multi-qubit systems with two-body interactions 

Let us now consider a system of N qubits with a Hamil- 
tonian H12...N made out of two-qubit interactions, 

H 12 ... N = J2 H v- ( 16 ) 

If interactions ify are arbitrary, the characterization of 
allowed simulations using H12...N and HLU becomes ex- 
cessively involved. However, in those systems where, e.g., 
Hij = H for some fixed interaction H , we automatically 
recover the Results 1 and 2 discussed above. This also 
happens if some couples of qubits do not interact at all. 

As an example, we may consider an array of spins 
with first-neighbour interaction o~ z ® a z . By applying 



A. Synthesis of symmetric two-qubit gates 

In what follows we take advantage of a result of |||o) , 
which says that any unitary matrix U acting on the 
Hilbert space C 2 ® C 2 can be written as 

U = ua®u b exp^— i^2 ^'k a k® cr/^j va®vb, (18) 

the decomposition being unique if we make the restriction 

tt/4 > Ai > > |A^|. (19) 

The recipe for getting A' is also given in ^0|. When 
imposing Eq. ( p"7| ) we get that the onc-qubit unitaries 
ua and va must be equal to «b and vb respectively, but 
since now we cannot modify the sign of the coefficients 
Afc , a unique decomposition for U arises only if we require 

7r/4>|Ai|>|A^|>|A^|. (20) 

Because we assume HLU can be performed arbitrarily 
fast and we are only concerned with simulating times, 
simulating U is equivalent to just simulating 

U^, = exp^-i^A' fc cr fe ® o-kj- (21) 

The construction U^, is called the canonical form of U. 
Note that all commutators ® o~j, o~k ® o~k] vanish and 
that exp(— iir o~k®o~k/2) = —icrk®o~k is a HLU. This im- 
plies that, for any vector n = (^1,7^2,^3) with integer 
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components n&, Up +jrf ,, 2 = Up U n ^/2, which is equiva- 
lent to Up up to HLU. Accordingly, we can associate to 
U any one of the Hamiltonians pl[ 



3/2 



3 

£ 

k=l 



(22) 



Now, if we are able to simulate the action of one of 
these Hamiltonians for a time t' = 1, we can imple- 
ment the gate U on our two-qubit system, since U and 



are HLU-equivalent. In addition, because Hab 



and Hp +2Lf , are both symmetric, Result 1 sets the sim- 
ulation condition: 



3 ft 



-< c A. 



(23) 



That is, even if Hab is not able to simulate Hp with 
HLU, it may still be able to simulate Hp +vf ., 2 for some 
vector n, and thus perform gate U. Notice, however, that 
in the case when A is proportional to (1,1,1) (isotropic 
interaction) there is no integer vector ft that makes the 
majorization possible unless A' is itself proportional to 
(1,1,1). In the next paragraph we explore for which As 
the condition ( |23] ) is satisfied whatever the value of A' is. 

The sum of the components of A, SfeAfc, can be pos- 
itive, negative or zero. If S^Afc is positive, by choosing 
n large enough, we can make A' -I- %(n,n,n) as close as 
we want to a vector proportional to (1,1,1), which is al- 
ways majorized by any cA not proportional to (1, 1, 1). 
The same happens if E/cAfc is negative, but in this case, 
n must be large and negative. If S/cAfc — 0, cA only ma- 
jorizes vectors having also the sum of their components 
equal to zero, and, not always exists an integer vector ft 
making £/c(A' fc + irnk/2) = 0. We conclude that, inter- 
acting hamiltonians with E^A^ different from zero and 
A not proportional to (1,1,1) can be used to synthesize 
any symmetric gate. 

Notice, from Eq. (^3|) that increasing the modulus of 
A' + 7m/ 2 by changing ft carries an increase of the time- 
overhead c (equivalently, decrease of the efficiency s of 
||). Thus, we can conlcude that the more isotropic Hab 
is (i.e. the closer A is to a vector proportional to (1, 1, 1)) 
the less efficient Hab is to synthesize gates. As a limiting 
case, if A is precisely proportional to (1, 1, 1), it will never 
majorize any A' not proportional to (1,1,1). 



B. Time-optimal gate synthesis 

So far we have seen that any symmetric gate can be 
synthesized by almost all interactions. A priori, the best 
synthesis protocol ought not to make use of a unique 
simulated interaction since there is freedom to change it 
along the protocol. Yet, the following result states that, 



indeed, time-optimal gate synthesis is achieved maintain- 
ing a unique simulated Hamiltonian. As before, A and A' 
are the vectors associated to Hab and U, respectively. 

Result 3 (Theorem 10 in |^|): The time-optimal 
way to synthesize a symmetric gate U using the inter- 
action Hab an d fast HL U consists of simulating, for a 
time t' — 1, among all Hamiltonians Hp,^- such that 



(A' + -ft) < cA 



(24) 



the one with smallest time overhead c. The minimal in- 
teraction time (interaction cost of JT^/j is given by c. 

The proof of this result follows from particularizing 
theorem 10 in Ref. || as indicated in p2|| , when taking 
into consideration that, up to permutations of the coef- 
ficients Afe |21|, Eq. (22) contains all Hamiltonians H' AB 
such that U is HLU-equivalent to exp(—iH' AB ). 

Let us illustrate the above construction with a spe- 
cific example. We would like to synthesize the non-local 
part of the CNOT gate, \i)A\j)B —* ®j)b, by us- 
ing a physical system which interacts according to the 
anisotropic exchange interaction: H = a x ®a x + o~ y ®o v , 
that is A = (1, 1,0). We first note that the CNOT gate 
is equal to: 



OE/Hadamarde* 
(/®C/Hadamard) 



that is, its interaction part reduces, up to HLU, to 



u 



(tt/4,0,0) = exp(— i—a x ®(j x ) 



(25) 



(26) 



A direct simulation of this Hamiltonian cannot be achi- 
eved since the condition (ir/4, 0, 0) -< c (1, 1, 0) cannot be 
fulfilled for any c (vector (1, 1, 0) is already more isotropic 
than (7r/4, 0,0)). Instead, we try to simulate the Hamil- 
tonian -ff(7r/4, tt/2. tt/2) 1 by which we would also achieve 
the goal according to (p2]). Now the simulation is pos- 
sible, (7r/4, 7r/2, 7r/2) -< c (1,1,0), with time overhead 
c = 5tt/8. It can be seen that this is the time-optimal 
way to do this. A possible mixing of permutations of the 
components of A = (1, 1, 0) leading to (7r/4, 7r/2, 7r/2) is 





(27) 



Note that to permute the x-z(y-z) axes we can perform 
a rotation of 7r/2 along the y(x) axis. Summing up, in 
order to synthesize the non-local part of the CNOT gate 
we have to divide the time t = 57r/8 into small equal 
intervals e and, in each one of those intervals, we have to 
perform the following sequence: 



iHlc 



1 e * 



°~b) 



e -iH§« 



(28) 
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The gate is synthesized once this protocol is repeated for 
a time t — 5tt/8. A more detailed explanation of the type 
of manipulations needed above is given in . 

IV. HAMILTONIAN SIMULATION AND GATE 
SYNTHESIS WITH AN ASYMMETRIC 
INTERACTION AND HLU 

Let us now consider how to simulate symmetric inter- 
actions and synthesize symmetric gates by using HLU 
and a two-qubit interaction Hab which is neither sym- 
metric nor antisymmetric under exchange of qubits. 

Recall that we can associate to Hab a real matrix M 
(its Pauli representation) through Eq. (|J). In what fol- 
lows we give a method which gets rid of the antisym- 
metric part of Hab, M a = {M — M T )/2, and the local 
terms, Ka®I and I®h,B, without affecting the symmetric 
part M s = (M + M T ) /2. To achieve this goal we use the 
following property: 

- o- k <g)a k (c„®cr TO ) u\®a\. = S nm a n (g)cr n , (29) 

where n and m take the values 0,1,2,3 and <To is the 
identity matrix. Performing this unitary mixing on a 
general Hamiltonian H, we get 

3 ^ 3 
^2 -g Ok®Ok H al®a\ = ^{M) kk a k ®o- k . (30) 

fe=0 fe=l 

That is, only the diagonal part of M survives. Suppose 
that Q is the 5*0(3) matrix that diagonalizes M s , and u is 
the unitary matrix corresponding to Q via the morphism 
(0). Then, if instead of (|3C|), we perform the mixing 

J2l=0 I {u ] <7 k u)®(u^(7 k u) H (u t <7fcu) t (g) (uVfcu)* 

= EL-=iW«^®oj. ( 31 ) 

we get the projection onto the symmetric interaction part 
of the Hamiltonian and no local terms, with no loss of 
efficiency In order to make a simulation using only the 
symmetric interaction part, we divide every infinitesimal 
time between the action of two HLUs into four equal 
intervals and perform the mixing @. Therefore, 

Result 4: Results 1 and 3 also hold in the case when 
the simulating Hamiltonian Hab is n °t symmetric. 



V. UNIVERSAL QUANTUM COMPUTATION 
WITH AN HYBRID MODEL OF SIMULATION 

We have already argued that an interaction and HLU 
are not sufficient to perform universal quantum compu- 
tation, since not all gates can be achieved using only 



these resources. In this last section we analize which ex- 
tra resources are needed so that the HLU model recovers 
universality. As extra resources we consider inhomoge- 
neous LU. A context where this hybrid model may apply 
is the case of two physical qubits which are brought very 
close to each other so that they can interact. When the 
qubits are close, only HLU can be performed in order to 
control the evolution, and the HLU model applies. But 
when the systems are brought far apart and no longer in- 
teract, they can be addressed individually and therefore 
arbitrary LU can be applied. 

We have seen in sect. Ill that any symmetric interac- 
tion (a part from the isotropic interaction and interac- 
tions with txM = 0) can be used to perform any sym- 
metric gate. This result has been extended in sect. IV 
also to asymmetric two-qubit interactions (with M s not 
proportional to the identity and trM ^ 0). Bearing this 
in mind, we announce the following result. 

Result 5: Any two-qubit gate can be performed using 
HLU and any given interaction Hab (such that M s is 
neither trace-less nor proportional to the identity), pro- 
vided that we allow for two non-homogeneous LU. 

To justify the last assertion, we recall that any two- 
qubit gate U can be written as in Eq. (|Tg|). Thus, we 
can use the available interaction Hab to perform the non- 
local part of U , which between two non-homogeneous LU 
interventions leads to U. 

This method is very inefficient when M s is almost 
proportional to the identity matrix or has very small 
trace, but alternative methods can be developed in these 
specific cases, and the hybrid model turns out to be 
fully universal. For instance, the isotropic interaction 
A = (1,1,1) can be used to perform the gate given by 
A' = 7r/8(l,l,l). Then, we can compose two of these 
gates with a local unitary to achieve a (tt/4, 0,0) gate, 
that is the non-local part of a CNOT gate, which to- 
gether with LU is known to be sufficient for universal 
quantum computation. 

VI. CONCLUSIONS 

We have studied the simulation of non-local Hamilto- 
nians and the synthesis of non-local gates when only HLU 
can be used to control the available interaction. This 
model is motivated by the difficulties that often arise in 
solid state physics, where individual systems can not be 
addressed independently because they are very close to 
each other. 

Although notoriously less powerful than inhomoge- 
neous manipulation, homogeneous local unitary trans- 
formations are sufficient to establish simulations within 
broad classes of Hamiltonians arranged by their sym- 
metry. In particular, HLU can be used to dispose of 
the anisotropic terms of any interaction, as required in 
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Ji2[ . Symmetric gate synthesis can be performed using 
symmetric interaction Hamiltonians and LHU. Moreover, 
these resources can be supplemented with just a few ex- 
tra non-homogeneous transformations to fully simulate 
any two-qubit unitary transformation, and thereby allow 
for universal quantum computation. 
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